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Abstract 



Particle creation at the expense of gravitational field might be sufficient to explain the cosmic evolution 
history, without the need of dark energy at all. This phenomena has been investigated in a recent work by 
Lima et-al [5] assuming particle creation at the cost of gravitational energy in the late Universe. However, the 
model does not satisfy the WMAP constraint on matter-radiation equality [12] . Here, we have suggested a 
model, in the same framework, which fits perfectly with SNIa data at low redshift as well as early Integrated 
S^ • ■ Sachs- Wolfe effect on matter-radiation equality determined by WMAP at high redshift. Such a model, requires 

pH ' the presence of nearly 26% primeval matter in the form of baryons and CDM. 

I < 

1 Introduction 

Recently released 7- year WMAP data [I], [2] has found no trace of deviation from the standard A CDM model. 
i— — i, But the problem in connection with the cosmological constant remains unresolved. The vacuum energy density, 
as calculated by the field theorists, is some 10 120 order of magnitude greater than the cosmological constant A 
required by the cosmologists, to explain late time cosmic acceleration, which is of the order of Hq , Hq being 
*y-\ • the present Hubble parameter. So far, many alternatives to the standard A CDM model have been proposed 
£C) , and explored and as a matter of fact, all of these models have been found suitable to explain late time cosmic 
0^ ' acceleration. The problems associated with these models are, - they can not be distinguished from the standard 
A CDM model at one hand, and most of them are not suitable to explain the early Universe, on the other. However, 
A CDM model requires 26% of matter in the form of pressureless dust, out of which only 4% are baryons and the 
rest, about 22% are cold dark matter (CDM). Since dark matter interacts only with the gravitational field, so it 
plays a key role in the structure formation. The gravitational Jeans instability allows compact structures to form 
and is not opposed by any force like radiation pressure in the case of dark matter. As a result, dark matter begins 
to collapse into a complex network of dark matter Halos, well before ordinary baryonic matter, which is impeded 
by pressure force. Without dark matter the epoch of Galaxy formation would have occurred at a substantially 
later stage, than observed. Thus, the amount of CDM (22%) must have been created in the very early Universe, 
prior to the radiation dominated era together with the baryons, by some sort of mechanism, viz., supersymmetry 
breaking, cosmic string decay or particle creation at the expense of gravitational field. These particles are usually 
supposed to be weakly interacting massive particles (WIMP). For example, as a heavy stable particle, the lightest 
neutralino is an excellent candidate to comprise the Universe's cold dark matter. In many models (see [3] for 
a nice review) the lightest neutralino can be produced in the hot early Universe and leave approximately the 
right relic abundance to account for the observed dark matter, ie., 22% as required by A CDM model. Now, if 
phcnomenologically one considers that CDM may also be produced by gravitational particle creation mechanism, 
even at a very slow rate, during the late time evolution of the Universe, viz., during the matter dominated era, 
as considered by Lima et-al [3], then it may be possible to explain the presently observable acceleration of the 
Universe, without taking dark energy into account. The very advantage of the creation of cold dark matter over 
dark energy is that it avoids coincidence problem and also may be detectable in future experiments. In the present 
work our focus is on the cosmological consequences of particle production on the evolution of late stage of the 
Universe which was initiated recently by Alcaniz and Lima [4 and Lima, Silva and Santos(lss) [5J. 

We remember that in the eighty's, the motivation that initiated to go after inflationary scenario, was the funda- 
mental three problems associated with Friedmann model, viz., the horizon, the flatness and the observed isotropy 
and homogeneity in the cosmic microwave background radiation (CMBR). There was an additional problem in the 
form of huge entropy per baryon in the observable Universe (~ 10 87 ) . Einstein's equations are purely adiabatic 



and reversible. Consequently, these equations can hardly provide, by themselves, an explanation relating to the 
origin of cosmological entropy. This problem was resolved [6] by taking into account the cosmological consequence 
of irreversible particle creation phenomena in the framework of Einstein's equation, classically. 

Particle creation phenomena was explored largely during the last century to explain the early Universe. Matter 
constituents may be produced quantum mechanically [7], [8], [9] in the framework of Einstein's equations. Cos- 
mological consequence of particle creation mechanism is studied taking into account an explicit phcnomenological 
balance law (see appendix) for the particle number [5], |10) . |11) in addition to the familiar Einstein's equations. 
In view of such a balance law, Prigogine et-al [6] successfully explained the cosmological evolution of the early 
Universe. 

Recently, Lima et-al [5] have developed a late time model Universe, taking into account the creation phenomena in 
the matter dominated era. The model admits early deceleration followed by a recent acceleration of the Universe 
as suggested by present observations and fits SNIa data to some extent. However, they [5] have shown that in 
their model the creation phenomena is never ending, as a result cosmic evolution does not ever include standard 
radiation dominated or matter dominated Friedmann era. This definitely creates problem in explaining the struc- 
ture formation of the Universe and the CMBR. Later, Steigman et-al [12] analyzed the model in the two limits of 
high and low redshifts. They have observed a clear conflict between the WMAP constraint on matter-radiation 
equality z eq at high redshift and SNIa data at low redshift. The main criticisms of the (3 — 7 model proposed by 
Lima et-al [5] are that they have not taken into account the amount of CDM created in the very early Universe at 
one end, and that their creation rate F = 3/3 H + 3"fHo depends on the present Hubble parameter, on the other. If 
these problems are alleviated, then phenomenologically particle creation process obviously unifies early inflation 
with late stage of cosmic acceleration in an elegant fashion. 

In the present work, we propose a model where, instead of choosing the creation parameter T arbitrarily, we 
have considered the experimentally verified fact that the Universe has recently entered an accelerated phase of 
expansion. As a result we have chosen the scale factor judiciously, such that particle creation could start again in 
the matter dominated era. This naturally alleviates the said problems and unifies the early and the late stages 
of cosmic evolution, in an elegant fashion. Additionally, the model fits perfectly with the WMAP constraint on 
matter-radiation equality z eq only if one considers the presence of nearly 26% of primeval matter in the form 
of baryons and CDM. In view of such a model, particle creation phenomena is now able to explain the history 
of cosmic evolution from the very early Universe till date, without requiring dark energy at any stage and thus 
avoiding coincidence problem. 

In the following section we write down the field equations incorporating the phenomena of particle creation, which 
is apparent through the balance equation. In section 3, we show how the conflict between high redshift and the 
low redshift data may be reduced in the model presented by Lima et al [5], just by accounting for some amount 
of CDM produced in the very early Universe. However, we also mention some more problems which are still 
associated with their [5] model. In section 4, instead of choosing a form of the creation rate arbitrarily, we rather 
choose a form of the scale factor, suitable for a transition from the early deceleration to the late time cosmic 
acceleration, to understand the associated problems. This gives us insight to find a form of the creation rate F 
associated with the so called intermediate inflation [13]. Such a form of T appears to be much elegant to study 
the cosmological evolution alleviating all the problems discussed. This has been done in section 5. Finally we end 
up with the conclusion in section 6 and an appendix in section (7), to calculate the balance law. 



2 Balance law and the Field equations 

As mentioned in the introduction, cosmological consequence taking into account particle creation phenomena is 
studied by using an explicit phenomenological balance law for the particle number [6], [10]. Such a balance law 
in the process of particle production is modeled by (nu a )- a = "J, with a source term ^ . For a vanishing vp the 
particle number is conserved. In the isotropic and homogeneous background metric, 
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above balance law reduces to, ^ = ^ + @ = ? = I\ where, u a , N , n and iJ are the fluid four velocity vector, the 
total number of particles, particle number density and the Hubble parameter respectively, while, = u a - a — 3H 
is the expansion scalar. Thus, the field equations in the spatially flat (k — 0) Robertson Walker metric (1) 
associated with particle creation phenomena are, 

2H + 3H 2 = -8TrG( Pm +p cm ), (2) 

3H 2 = 87rG( Pm + p cm ) = 8nGp, (3) 

3H+- = T=-, (4) 

n n 

P + Prn ,_. 

Pcm = — ^— r. (5) 

For the last equation (5) please see the appendix. In the above set of equations, H = - is the Hubble parameter, 
p m and p m are the pressure and the energy density of the matter existing in the Universe in the form of a 
barotropic fluid containing baryons and cold dark matter, created in the very early Universe. p cm and p cm are 
the pressure and the energy density of the cold dark matter in the form of WIMP created at the late stage of 
cosmic evolution, ie., during matter dominated era. T is the creation parameter, and n is the particle number 
density. It has been shown in the appendix that the second law of Thermodynamics allows the creation of particle 
from the gravitational field and the process is irreversible. Thus the creation parameter T > 0, and so it is clear 
from equation (5) that the particle creation phenomena is always associated with a negative pressure p C m , which 
may be responsible for acceleration at the late stage of cosmic evolution. We would like to mention at this stage 
that, other than cold dark matter and baryons, different types of particles may be created in view of quantum field 
theory in curved space time [7], [8], [9] and all are associated with a negative pressure, as mentioned. However, 
also, as mentioned in the introduction that the structure formation requires 22% that are constituted by dark 
matter, which must have been created in the early Universe. The present estimated amount of baryons is 4% and 
the rest amount required for present cosmic acceleration is usually treated as dark energy. Other form of matter 
(like hot dark matter, say) has negligible contribution. Here, we proceed to show that creation of the same amount 
of CDM also solves the puzzle. Further, note that the creation of baryons and CDM in the early Universe was also 
associated with a large negative pressure. However, the creation phenomena weakened and finally stopped, when 
Universe expanded sufficiently, thereby giving way to the hot big bang followed by the radiation dominated era of 
Friedmann type (a oc t* ) [6]. There after, the baryons and the CDM created in the very early Universe obviously 
start acting as pressureless dust, so that we can take p m ~ and hence p m = poa~~ 3 , po being a constant. Thus 
we can simplify the above equations to get, 

r = 3H + 2 | = ZH ( 2H ^ H l H2 ) = -Mw e , (6) 

3-ff _ 3 . . 

Pcm = g-g - Poa - , (7) 

Pcm = ~^c HT > ( 8 ) 

w e being the effective state parameter. In view of the above set of three equations (6) - (8), we need to find 
the scale factor (a) (and consequently (H) , the Hubble parameter), the creation rate T, the creation pressure 
p cm and the creation matter density p cm . Obviously, we need yet another suitable condition to solve the system 
of equations. Lima et al [5] studied these equations under the assumption of a form of the creation rate T . 
Afterwards, while further studying their model in connection with data fitting, they found a clear conflict between 
the low (SNIa) and high (WMAP constraint on z eq ) [12] redshift limits. In the following section, we review the 
problem and show how the production of CDM in the very early Universe alleviates the problem. 



3 A brief review of lss model 

The Friedmann equation, taking into account the created matter, baryonic matter and radiation reads, 
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where, p c is the present value of critical density. To calculate the last term let us take the total number of created 
particles at an instant to be N = nV , where V = Vo(l + z)~ 3 is the comoving volume at that instant. Thus the 
creation rate is given as, 
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Jt 
where, p cm o is the present value of the created matter density. Thus, 

^m = n cm (i + zf ex P (- / ° r*'), (12) 

Pc Jt 

where, f2 cm is the density parameter corresponding to the created matter. So, the Friedmann equation (9) finally 
reads, 

H \ 2 /.to 



— =fi r (l + z) 4 + f7 B (l + z) 3 + f7 cm (l + z) 3 exp(- / rdt'). (13) 

Now, under the assumption T = 3/3 iJ + 3jHq , where, /3 and 7 are constants and Hq is the present Hubble 
parameter, (lss) [S] obtained a solution of the scale factor in the form, 



a(t) = a 
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which admits the observed transition from early deceleration to late time acceleration. In a later investigation |12j , 
this model was found to produce a conflict between SNIa data at low redshift and WMAP - 5 year data constraint 
[M] on mater-radiation equality z eq — 3141 ± 157, occurred at the high redshift limit of observed ISW effect. Let 
us review the situation to find the real problem associated with the conflict. The Friedmann equation (13) in the 
model under consideration reads, 



Upj = n r (l + zf + n B (l + zf + fl cm (l + z)^ 1 '^ exp3 7 (r - r ), 



(14) 



where, r = H$t and tq — iJo^o are the age at any instant and the present age of the Universe respectively, in the 
units of Hubble's age (Hq 1 ) ■ Setting, Q cm = 1 — Qb , the 7 — f3 relation is obtained (see equation (34) in [12]) as, 

7 = (1 - /?) [(1 - n B f - {a r (i + Zeq ) - n B }Hi + z eq ) 3 ^ . (15) 

This model fits SNIa data for /3 = and 7 = 0.66 ± 0.04, while 1 + z eq = 1798±f|| , taking fi B = .042. Clearly, 
the model does not fit with the WMAP - 5 year data constraint [TJ on mater-radiation equality z eq = 3141 ± 157, 
occurred at the high redshift limit of observed ISW effect. This contradiction may be alleviated easily, if we consider 



existence of cold dark matter that was created in the early Universe and which was responsible for inflation. As, 
already mentioned, this amount of CDM created in the very early Universe behaves now as pressureless dust and 
has been redshifted like baryons. If we now add corresponding density parameter CIqdm , associated with the cold 
dark matter created in the very early Universe in equation (13), it reads 

JH =ft r (l + z) 4 + f! m (l + 2 ) 3 + fl cm (l + z) 3 exp(- J°Tdt'), (16) 

where, Q, m = fls + Qcdm and, £7 cm = 1 — Q m . In the absence of matter creation phenomena in the late Universe, 
/3, 7 vanish, and hence T = 0. Thus there is no creation pressure p m as well as creation matter density p cm - 
Hence, Q cm = 0. Thus at the matter-radiation equality (z eq — Kj®- — 1) taking, Q m = fls + &cdm — 0.26 and 
fl r = 8 x 10~ 5 , one recovers z eq = 3249, which is at par with WMAP data. Equation (15) now takes the form, 

J = (1-/3) [(1 - m )5 - {Q r (l + Zeq ) - fl m }?{l + z eq )^ . (17) 

If we now consider that 16% of CDM (say) were produced in the very early Universe, then fl m = ^Ib+^cdai = 0.2 
and thus for (i = and 7 = 0.66 ± 0.04, z eq = 3186±2i4, which is very much at par with WMAP data p], [2], 
|14j . This clearly indicates that one should include the contribution of CDM created at the very early Universe. 
The creation of CDM in the very early Universe, as mentioned in the introduction, was halted and the Universe 
entered usual Friedmann radiation dominated era. Thereafter, this amount of CDM is being redshifted like baryons. 

In the above analysis, while we have showed how the conflict encountered between low and high redshift data |12j 
may be reduced, nevertheless, it does not support the model [5]. Firstly, in their (3 — 7 model, /3 = 0, somehow 
fits SNIa data, which is not a very good fit at all (see fig. 1 [H]). Further, f3 — turns out to give a constant 
creation rate throughout the evolution of the Universe, which is highly objectionable. Also, the choice of the 
creation parameter ( T ) as a function of present Hubble parameter ( H q ) implies that the model is plagued by the 
coincidence problem. Finally, we could accommodate only 16% of CDM out of 22% to alleviate the conflict [12) . 
Addition of another 6% of CDM shifts z eq to a much higher value. In view of the above criticism we pose to 
present a more realistic model. 

4 Case 1 

To get an explicit solution of the field equations (6) through (8), we need yet another physically reasonable 
assumption. However, we really have no idea of the rate of matter creation T cither from quantum field theoretic 
(QFT in CST) or from classical kinetic approach. Nevertheless, it is clear from equation (5) that such phenomena 
is associated with a negative pressure p c . If the creation pressure is sufficiently negative, it might lead to an 
accelerating phase of cosmological evolution. Additionally, in order to get an idea about the form of the creation 
rate we can also depend on the presently available cosmological data and the best fit models. Since A CDM 
model has excellent fit with the SNIa and WMAP data, so we can infer certain important aspects of cosmological 
evolution. First of these is definitely that the Universe has encountered a transition from early deceleration to late 
time acceleration. Next is, early growth of perturbation should track A CDM model closely. These facts allow us 
to choose a suitable form of the scale factor at par with the present experimental results and in the process, we 
expect to get an idea on the form of T . In view of the two aspects of the late time cosmological evolution just 
discussed, we choose the scale factor in the matter dominated era as, 

a = a ti +b t a =b {qt^ +t a ), (18) 

where, ao, b and a > 1 are constants with, q = #£ > 0. With such a form of the scale factor, the first term 
of equation (18) plays the leading role in the early stage of cosmological evolution in the matter dominated era 
and the Universe tracks Friedmann model a oc t$ , closely. Thus, the growth of perturbation in connection with 
the structure formation tracks A CDM model closely. The second term appears in connection with the creation of 
matter, which is associated with a negative creation pressure. As t increases, second term starts playing significant 
role and it starts dominating as t a ~~z > q, which implies that the creation phenomena starts rather late. Finally, 



since a = bo[ r + oc(a — l)t a 2 ] , so acceleration (a > 0) starts, only when t > 
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Thus, equation 



(18) clearly depicts early deceleration and late time acceleration of the Universe. Now, in the matter dominated 
era, p m = , the solutions in view of the chosen form of the scale factor (18) are, 
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(24) 



qti + t a 

In the above set of solutions w e and w cm are the effective state parameter and the state parameter corresponding 
to the created matter respectively. Since a > 1 , so the expression (23) for w e is clearly negative which, as men- 
tioned earlier, implies that the created matter is such that it is always associated with a negative creation pressure. 



Fitting the observational data. 

The present model is parametrized by the two parameters a and q. With a = 4, h = -|p|- = 0.7 Gyr^ 1 , and 

fixing Hoto = 1 , q is automatically fixed and the distance modulus versus redshift curve (blue) is found to fit 
perfectly with the ACDM model (red). In fact the two are practically indistinguishable (fig.l). We get the present 
value of the effective state parameter w e o = — 1 and the transition redshift z a = 0.56, which are in excellent 
agreement with ACDM model. It is observed that the effective state parameter encounters a transient double 
crossing of the phantom divide line in the future (fig. 2) and so Big-Rip singularity is bypassed. For a < 3.2, the 
phantom divide line is never crossed, while for a > 4 , the first crossing occurs in the past but second one always 
occurs in the future. The model fits perfectly with SNIa data for a wide range of values of a . One can also 
observe that the state parameter w e remains nearly zero till z = 2.5 (fig-2), which confirms that the growth of 
perturbation in the present model tracks the concordance model closely. 



Thus the present model fits perfectly with the ACDM model without any problem what-so-ever. However, the 
problem arises while one tries to fit the recently released 7-year WMAP [1] , [2] constraint on the redshift of matter- 
radiation equality at early ISW effect. The value of the integral in equation (16) is X = exp (— J ° Tdt r ) = 0.3236 , 
taking a = 4 . As a result the redshift at matter-radiation equality is pushed far away to z eq = 4380 , taking only 
baryons into account, ie., fl m — Qb = 0.042. It goes even further if some amount of Qcdm is incorporated. The 
situation is even worse for a > 4 . The whole situation taking some lower values of a and with the same values of 
h = 0.7, H t — 1 , is depicted in table - 1. 
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Figure 1: Distance modulus (M — m) versus redshift z plot of the present model (blue), shows perfect fit with 
the ACDM model (red) for a = 4. In fact, for a > 1.5 the fit is perfect and it remains so, upto a = 200. 




Figure 2: State parameters w e (z) (continuous) and w cm (z) (dashed) have been plotted against the red-shift 
parameter z for a — A. Both remain nearly zero till z — 2.5. The transition redshift is z a rs 0.56. The present 
value of w e is — 1. A smooth transient crossing for w cm is observed. The first crossing occurs in the past while 
the other will occur in the future. For a < 3.2 , no such crossing is observed. 



Table - 1 



a 


X 


Nature of T 


Za 


W e0 


J l ni 


&CDM 


WcmO 


%eq 


Fit with SNIa 


4 


0.3236 


Rises from z — 30 


0.56 


-1 


0.04 


Nil 


-1.20 


4380 


Indistinguishable from ACDM. 


2 


0.1876 


Rises from z — 2500 


1.08 


-0.56 


0.09 


0.05 


-0.62 


3257 


Indistinguishable from ACDM. 


1.8 


0.1563 


Rises from z = 7000 


1.2 


-0.51 


0.12 


0.08 


-0.58 


3217 


Indistinguishable from ACDM. 


1.5 


0.0962 


Very large initially 


1.46 


-0.44 


0.18 


0.14 


-0.545 


3234 


Fit is not the very best. 



We are in search of a model where additional creation starts some time after the matter-radiation equality z eq . 
r — z plot (figure 3) depicts that the creation rate is almost vanishing (~ 10 -8 ) till z = 30 for a — 4 and only 
thereafter it rises sharply. However the problem is that the creation rate is not sufficiently large to fit WMAP 
constraint on z eq . For a > 4, the situation is even worse. The behaviour is the same also for lower values up-to 
a = 2. The only difference is that the creation starts earlier in this case at z = 2500, and so created matter is 
a little large. Thus WMAP constraint on z eq is satisfied taking into account only a small amount of CDM, viz., 
^cdm — 0.05 . Thus the problem with lss model is encountered here too. For a < 2 , creation starts very early, 
much before z„„ and so these cases are discarded. 
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Figure 3: Creation rate T versus red-shift parameter z for a = 4, shows that it is nearly vanishing in the past and 
started developing only recently at z « 30 . For higher values of a , T starts developing from even smaller redshift 
value, while for a < 4, it starts earlier. Thus, WMAP constraint on z eq can be fitted for a < 3.5, accommodating 
some amount of CDM. 

5 Case 2 

The main problem, we repeat, associated with the above model and the lss model [5] is that the creation of 
matter is not sufficient to fit WMAP constraint on the redshift of matter-radiation equality at early ISW effect, 
accommodating 22% that are constituted by dark matter created in the very early Universe. Thus, the choice 
of r should be such that the creation of matter starts near z « 3000 and in the later epoch it should increase 
considerably, to make X sufficiently small. So, to find a suitable form of T , we try with a scale factor associated 
with the so called intermediate inflationary solution [13 , viz., 



ao exp 



[At% 



(25) 



a being a constant. Such a solution for A > and < / < 1 , was presented by Barrow[T2], and was shown to 
lead to late time acceleration [TS] in different models. To appreciate the underlying beauty of the ansatz (25), let 
us expand it as, 



a = a [l + At f + —A 2 t 2f + ], 

and observe that for / = § , the standard matter dominated era of Friedmann model is recovered in the early 
Universe when second term dominates, and the third term becomes responsible for accelerated expansion in the 
late stage of cosmological evolution. For / = | , the third term leads to the standard Friedmann model and 
acceleration starts a little late. For even smaller values of /, the model tracks decelerated expansion for a longer 
time recovering the standard Friedmann model at some intermediate stage of evolution and leads to accelerated 
expansion at much later stage of cosmic evolution. 
The redshift parameter 1 + z = ^jfr , where, to is the present time, is found as 



l + z = cxp[A(^-i / )]. 
Hence, the Hubble parameter takes the following form, 
Af Af 



H = 
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(26) 
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where, we have used equation (26) to get the second equality in the expression of H . In view of equation (6), we 
can now find a form of T as, 



r = 3H - 2(1 - /) 



H 
Af 



(28) 



This form of T is clearly different from the /3 — 7 model [5] . The most important difference is that the creation 
rate T here, starts developing only when the Hubble parameter, 



H > 



2(1-/) 



(1-/) 



Af 



(29) 



since, as already mentioned in the introduction, T < is not allowed by the second law of Thermodynamics (see 
appendix). The creation pressure and the creation matter density are now found as, 



8irG Pcm = -TH = H 
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PmO 



3H 2 (l + z) 3 = 3H 2 n m (l + z) 3 , 



(32) 



in which, p c , p m o and ft m are the present values of critical density, matter density and the matter density 
parameter respectively. We can find the effective state parameter and also the state parameter of the created 
matter as, 



2H + 3H< 
ZH 2 



= - i+ Kw) 



(33) 
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(34) 



Now let us see how far this model parametrized by the two parameters A and / fits with the observed data. 

We have kept 0.96 < H Q t < 1 and 0.67 < h{= $M Gyr- 1 ) < 0.7, as par with HST project [TB]. As already 

-"0 
mentioned, the restriction on the parameters are, A > and < / < 1 . We have tested the model by choosing 

A and / which fit SNIa data, from a wide range of values between 0.08 < A < 25 and 0.03 < / < 0.99. The 

fit requires large A for small / and vise- versa. We have presented our results briefly in the following table - 2, 

taking only some integral values of A starting form A = 15 , since for lower values this model does not probe to 

large redshift z. We have taken z eq = 3300, which is very much at par with recently released WMAP data Q], 

[2] and VLb — 0.042, to find the amount of matter produced in the late stage of cosmic evolution restricting the 

amount CDM produced in the very early Universe. 



Fitting the observational data. 



Table - 2 



A 


/ 


zr=o 


Za 


WeO 


i L r!l 


i L crn 


QcDM 


15 


0.056 


468 


0.70 


-0.35 


0.246 


0.754 


0.204 


16 


0.053 


658 


0.71 


-0.35 


0.249 


0.751 


0.207 


17 


0.051 


1145 


1.32 


-0.36 


0.255 


0.745 


0.213 


18 


0.048 


1350 


0.82 


-0.35 


0.255 


0.745 


0.213 


19 


0.046 


2050 


1.04 


-0.36 


0.257 


0.743 


0.215 


20 


0.044 


2914 


1.09 


-0.36 


0.258 


0.742 


0.216 


21 


0.0415 


3080 


0.40 


-0.34 


0.257 


0.743 


0.215 


22 


0.0392 


3157 


-0.11 


-0.33 


0.256 


0.744 


0.214 



In table 2, zr=o an d z a symbolize the redshift values at which the creation of matter and the acceleration start 
respectively, while w e o is the present value of effective state parameter. Let us list our observation point by point. 

1. The distance modulus versus redshift curve fits between the present and the ACDM model (taking fl\ = 0.74 
and f2 m = 0.26) almost perfectly for a wide range of values of the parameters A and /. 

2. It is observed that for the combinations of A and /, which can probe to a distant redshift, the present value of 
the state parameter is nowhere near —1, yet, the model fits both the experimental data, viz., SNIa and WMAP. 
Particularly, for A > 22 , the acceleration is yet to start. 

3. The most important point is to note that for A > 10, z eq is at par with the recent 7- year WMAP data pQ, 
[2], only if 24% — 26% of matter (baryons and CDM) is assumed to have formed in the very early Universe. The 
table shows that the density parameter 0.74 < £l cm < 0.76, which corresponds to 74% — 76% of matter created 
in the matter dominated era. Thus, instead of taking into account 74% of dark energy, creation of dark matter 
by the same amount in the matter dominated era, solves the cosmic puzzle. 

4. The behaviour of the creation parameter T given in equation (28) has been plotted in figure (4) for a particular 
pair of the parameters A = 19 and / = 0.046. It shows that the creation started at z = 2050, reaches a maxima 
during reionization era and presently it is insignificantly small. The behaviour is the same for all other pairs of A 
and /, only the redshift values at which creation starts (z-p — 0) and its maxima changes. 

5. Taking the same values of A and /, figure (5) has been plotted. It represents the combined plot of effective 
state parameter w e and the state parameter w cm corresponding to created matter, versus redshift parameter z , 
since creation started. These two figures (4) and (5) depict that though creation started rather early at z — 2050 
and reaches its maxima around z — 1100, acceleration started only recently at z = 1.04. Using relation (26) it 
is found that it requires nearly 13.99 Gyr to create 74% of matter. On the other hand, inflation is supposed to 
start at 10~ 42 s. , and ends at around 10 _32 s. Thus 22% that are constituted by dark matter has been created in 
10 _32 s only, in the very early Universe. This gives a comparison of creation phenomena in curvature dominated 
and low curvature regions. 

6. We have also presented a suitable contour plot in figure (6), to explore the data presented in table-2 at a glance. 
The plot presents all the successful combinations of the parameters A and / , which fit SNIa data, and satisfy 
WMAP constraint on z eq = 3300, keeping H t « 1, and 24% < fl m < 26%. Calculation shows that WMAP 
constraint on z eq is not satisfied for lower or higher values of fl m , with the same parametric combination of A 
and /. Particularly, for Vt m = 0.2, 2400 < z eq < 2500, while for Q m = 0.3, 3800 < z eq < 3900. Thus, nearly 
26% of primeval matter in the form of baryons and CDM, is required to fit presently observable data, in view of 
particle creation phenomena. 

Finally, it is no less important to understand if adiabatic process occurs instead, due to the fact that large time 
taken to (eg., about 14 Gyr.) particle creation phenomena in the low curvature region. It is known [8] that if the 
expansion rate is very weak, the production of high mass particles is exponentially small. This is due to the fact 
that large amount of energy must emerge from changing gravitational field to supply particle's rest mass. Thus 
particle number remains adiabatic invariant and a comoving particle detector remains unexcited, which means 
the probability of detecting particles falls sharply to zero. However, in discussing quantum particle production in 
curved space time [8], the balance law (5) has never been accounted for, which is crucial in the present analysis. 
The balance law introduces a back reaction phenomena. As soon as some particles are produced, they impart 
negative pressure p cm , which enhances the expansion rate causing more particles of higher mass to produce. The 
process continues as long as the Universe expands sufficiently so that the curvature fluctuation is further reduced 
and the creation rate falls. This feature is present in the T — z plot of figure (4). Thus, it appears that adiabatic 
process will not be simulated. 



6 Concluding remarks 

We have presented a phenomenological cosmological model based on particle creation in the matter-dominated era, 
which fits SNIa data and the redshift of the early integrated Sachs- Wolfe effect at the matter-radiation equality, 
z eq = 3145±^ [I] and z eq = 3196i^g [2], determined by WMAP. The value of h-parameter (h « 0.70) and 
iJoio ~ 1 also are very much within observational limit obtained from HST project [16| . Such a model constraints 
the amount of primeval matter to 0.24 < Q m < 0.26, out of which CDM amounts to 0.20 < flc'DM < 0.22, which 
is again the same amount required for structure formation. 
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Figure 4: The behaviour of V versus z has been depicted for A — 19 and / = 0.046. Creation starts in the matter 
dominated era around z = 2050 and its rate has a maxima around zr mam = 1100(= z rec ombination) ■ Presently the 
creation rate is insignificantly small. The behavior is the same for all other combinations of A and / , only zr=o 
and zr max are different. 
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Figure 5: The figure shows the combined plot of effective state parameter w e and the state parameter w cm of 
created matter, versus redshift parameter z, since creation started, taking A = 19 and 0.046. While figure-4 
depicts that creation started at z = 2050 and reaches its maxima at z = 1100, figure-5 shows that most of the 
time universe undergoes decelerated expansion while acceleration started recently at z = 1.04. 




Figure 6: The contour plots of Hoto — 1 (dotted line) and il m for different parametric values of A and / 
which fit SNIa data, have been combined together. The plot shows that Hgto = 1 line lies within the two lines 
Q. m = 0.24 and fl m — 0.26, which are calculated taking z eq = 3300. Thus, the present model fits SNIa data, 
satisfies WMAP constraint on z eq , demanding the correct amount of il m required for structure formation and 
agrees with experimental constraint on H$to = 1 . 



Quantum particle production phenomena (QFT in CST) has been discussed in the literature in detail [8], [9]. 
The energy of these particles may then be extracted from the gravitational field [17]. To study the classical 
consequence of particle creation phenomena, kinetic collision theory may be adopted to find a balance law in 
addition to standard Einstein's equations, if weakly interacting massive particles (WIMP) are taken into account. 
In view of such a balance law a cosmological model of the early Universe has been explored [iaj . In that model, as 
the particle production rate becomes comparable to the expansion scalar ( — 3H ) , it builds up a large negative 
creation pressure that pushes the Hubble parameter H to approximately constant value. As a result, inflationary 
behaviour due to a large particle production rate is realized. Consequently, the universe starts with a de Sitter 
phase avoiding cosmological singularity. As inflation continues, the expansion rate becomes too large comparable 
to the particle production rate. In such a dilute cosmic fluid, particle production rate decreases and halts at some 
time. At this stage, in the absence of sufficient negative pressure, inflation ends giving way to reheating and the 
resulting Universe smoothly approaches the familiar Friedmann-Lemaitre-Robertson- Walker behaviour. Since the 
radiation dominated era, in such a model is of standard Friedmann type (i.e., a oc t 1 ' 2 , a being the scale factor), 
so the standard Big-Bang-Nucleosynthesis (BBN) remains unaltered. Thus cosmological evolution of the early 
Universe may be explained successfully in view of particle creation phenomena. 

Inflation makes the Universe almost spatially flat, which has also been confirmed by recent observations. However, 
at the end of inflation, if the Universe remains slightly away from spatial flatness, it may cause particle produc- 
tion in the matter dominated era again. If this happens, then even a slow particle production rate may cause 
sufficient negative pressure in billions of years to cause recently observed cosmic acceleration. This phenomena 
has been studied earlier 0], [5], but the model suffered from a clear conflict [H] between the low (SNIa) and 
high redshift (WMAP) data. This problem has been resolved in the present model, which constraints primeval 
matter to O m ~ 26% . Since, particle production starts long after matter-radiation equality, so the early growth of 
perturbation in connection with the structure formation also remains unaltered following the A CDM concordance 
model closely. Thus particle production process may successfully explain the late stage of cosmic evolution also. 
As a result, the cosmological evolution of the Universe may be explained successfully in view of particle creation 
phenomena from early time till date. 

Presently we are having numerous dark energy models, explaining the late time cosmic phenomena. Most of these 
models do not explain the early Universe at one hand and it is not possible to identify these models from one 
another in any of the future experiments, on the other. Cosmological consequence of particle creation phenomena 
does not require dark energy at all at any stage, and some programmes have been taken in the recent years to 
detect lightest neutralino of roughly 10 — 10000 GcV - one of the weakly interacting massive particles (WIMP). 
Thus, if creation phenomena of cold-dark-matter can solve the presently observed cosmic puzzle single handedly, 
without taking into account the dark energy at all, may be resolved in near future. 



7 Appendix: Thermodynamics of adiabatic particle creation 

In this appendix we formulate the balance equation in connection with particle creation phenomena. This has been 
done in [5], [BJ. However, the approaches are slightly different, hence we produce a straight forward calculation. 
Adiabatic cosmological evolution in the presence of particle creation can be treated in the open system, and so 
the first law of Thermodynamics is modified as, 

d{pV) + Pm dV - -d{nV) = 0, (35) 

n 

where, p, p m , V, n and h are the total energy density, the true thermodynamical pressure, any arbitrary co- 
moving volume, the number of particles per unit volume and the enthalpy per unit volume respectively. In the 
case under consideration, the system receives heat only due to the transfer of energy from gravitation to matter. 
So, creation of particles acts as a source of internal energy. Thus for adiabatic transformation the second law of 
Thermodynamics reads, 

TdS = d(pV)+ Pm dV - nd{nV) , (36) 

Combination of the two laws (35) and (36) gives, 

TdS = -d(nV) - ndinV) = TcrdN, (37) 

n 

where, we have used the usual expression for the chemical potential as pn = h — Ts and define s = £ to 
be the entropy per unit volume and a = ^ as the specific entropy. Thus we observe that the second law of 
thermodynamics viz., dS > implies dN > 0, and the reverse process is thermodynamically impossible, ie., 
particle can only be created and can not be destroyed. Further, expressing S in terms of a , the above equation 
can also be expressed as, 

TNda = => a = 0, (38) 

Hence, in the adiabatic particle creation phenomena, entropy increases, while the specific entropy remains constant. 
First law given by equation (35) can also be expressed as, 

Vdp + pdV + Pm dV - hdV - —dn = => Vdp - —dn = => p = h-, (39) 

n n n 

Now, the energy-momentum tensor T^ along with the conservation law when creation phenomena is incorporated 
are, 

T^ - {p + p m + Pcm )u^u v - (pm+PcmW, T£ tlt - 0, (40) 

where, p — p m + p cm is the total energy density and p m is the thermodynamic pressure, as already stated, while, 
p cm is the creation pressure and w M is the component of four velocity vector. The energy conservation law (40) 
in homogeneous cosmological models reads, 

P + @(p + Pm+Pcm)=0, (41) 

where, = ZH is the expansion scalar, H being the Hubble parameter. If we now plug in p from equation (39) 
in the above equation (40), we get, 

P + Prn / Q , h\ P + Prn ,,„, 

Pom- q- © + - = Q-r, (42) 



where, r = + - is the creation rate. 
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